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ABSTRACT 

We present an analytical model for the cosmological evolution of the FRII source 
population. Based on an earlier model for the intrinsic radio luminosity - linear size 
evolution of these objects, we construct theoretical source samples. The source dis- 
tributions in the radio power - linear size plane of these samples are then compared 
with that of an observed flux-limited sample. We find that the source parameters de- 
termining the radio luminosity of FRII objects can not be independent of each other. 
The best-fitting models predict the jet power to be correlated cither with the life time 
of the source or with the shape of the density distribution of the source environment. 
The latter case is consistent with the observed tendency of the most luminous radio 
sources at high redshift to be located in richer and more extended environments than 
their low redshift counterparts. We also find evidence for a class of FRII sources dis- 
tinctly different from the main population. These sources are extremely old and/or 
are located in very underdense environments. The luminosity function of FRII sources 
resulting from the model is in good agreement with previous results for high lumi- 
nosity sources. The apparent luminosity evolution of the radio luminosity function is 
not reproduced because of the high flux limit of the used comparison sample. The 
cosmological evolution of the median linear size of FRII sources is found to be mild. 

Key words: galaxies: active - galaxies: evolution - galaxies: jets - galaxies: lumi- 
nosity function, mass function - cosmology: theory 



1 INTRODUCTION 

This is the third in a series of papers presenting an an- 
alytical model for the evolution of FRII radio sources. In 
Kaiser & Alexander (1997), hereafter KA, we developed a 
model for the general dynamics of these sources. The intrin- 
sic evolution of the radio luminosity of a given radio source 
as a function of its linear size was investigated in Kaiser, 
Dennett-Thorpe & Alexander (1997), hereafter KDA. In this 
paper we use the model for the intrinsic radio luminosity- 
linear size evolution of FRII sources of KDA to constrain 
the cosmological evolution of these sources and their pro- 
genitors. 

The radio power (P) - linear size (D) diagram was in- 
troduced by Shklovskii (1963) as a powerful tool to inves- 
tigate the evolution of extragalactic radio sources. Baldwin 
(1982) pointed out that the P-D diagram is analogous to 
the Hertzsprung-Russel diagram for stars. However, there is 
a strong correlation between the radio luminosity and the 
redshift of sources in flux- limited samples. This Malmquist 
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bias implies that the source distribution in the P-D plane 
of such samples is a result of the intrinsic evolution of indi- 
vidual sources and the cosmological evolution of the source 
population as a whole. KDA point out that the scarcity of 
radio sources with linear sizes greater than roughly 1 Mpc 
can be partly explained by the steepening of the evolution- 
ary tracks through the P-D diagram of FRII sources caused 
by the energy losses of the relativistic electrons in their co- 
coons due to inverse Compton scattering of the Cosmic Mi- 
crowave Background Radiation (CMBR). Since the energy 
density of the CMBR increases with redshift, the cut-off in 
linear size above which only very few sources are found will 
put important constraints on the cosmological evolution of 
the radio source population. 

It is not clear why AGN formation occurs and how jet 
activity in galaxies is subsequently triggered. It is therefore 
also unknown whether all massive galaxies in the universe 
go through a phase of jet activity as part of their evolution 
or whether only those galaxies involved in violent processes 
like galaxy mergers can become active. 

For the purpose of the analysis presented in this paper 
we will simply assume that there is a population of pro- 
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genitors of radio sources without specifying exactly what 
they are. The only property these progenitors must have, is 
that they become active with a certain probability at some 
cosmological epoch, as described by the 'birth function', and 
then turn into a radio source whose intrinsic luminosity evo- 
lution is determined by the properties of the jets and the en- 
vironment the progenitors are located in. With this assump- 
tion there are then two possibilities how cosmic evolution in 
the number density of radio sources can occur. Either all pro- 
genitors existing at a given redshift become active at some 
point in their life time and the evolution reflects the num- 
ber of progenitors created at or before this redshift; or the 
progenitors are created very early in the evolution of the uni- 
verse, have very long life times and jet activity is triggered 
by a process not intrinsic to the progenitor. In this case any 
evolution in the number density of radio sources indicates 
that the probability for this process to occur is varying with 
redshift. For convenience we will use only the latter inter- 
pretation in the formulation of the models in the following 
sections. However, it should be born in mind that there is 
no way inherent in the model to decide which interpretation 
is the more likely. 

Assuming a birth function and a set of distribution func- 
tions characteristic of the environments of the progenitors, 
it is possible with the aid of the model described by KDA 
to predict a P-D diagram which accounts for observational 
selection effects. To constrain the models the predicted P-D 
diagram is then compared directly with the observed dis- 
tribution of sources. Knowledge of the selection effects, and 
hence the use of a complete observational sample is essen- 
tial for a proper statistical comparison. The only published 
sample of extragalactic radio sources with the completeness 
required for this analysis is the sample of Laing et al. (1983), 
hereafter LRL. In the original form the authors showed that 
the LRL sample is 96% complete to a flux limit of 10 Jy at 
178 MHz for sources with angular sizes less than 10'. Riley 

(1989) showed that the sample is complete even for extended 
sources down to the surface brightness limit of the 6C sur- 
vey of 120 mjy per beam (Baldwin et al. 1985). The sample 
subtends a solid angle of roughly 4.1 sr and contains 173 
radio sources of which 30 are of type FRI. The remaining 
143 sources are either quasars (43 sources) or radio galaxies 
(100 sources) with FRII morphology. We will assume here 
that quasars and FRII radio galaxies form one class of ob- 
jects distinguished only by a different viewing angle at which 
they are observed (e.g. Barthel 1989). The model used here 
for the calculation of the radio emission of FRII sources does 
not include the emission of the radio core. The low selection 
frequency of LRL ensures that relativistically beamed radio 
emission from the radio core of FRII sources does not con- 
tribute significantly to the overall flux. However, two sources 
in the sample, 3C 345 and 3C 454.3, may have been included 
in the sample only because of their boosted core emission. 
The changes to the statistical properties of the LRL sam- 
ple introduced by excluding these two sources are negligible. 
The source distribution of the LRL sample in the P-D plane 
is shown in Figure |l| 

Once the most likely birth function and distribution 
functions of the progenitor properties is found, it is straight- 
forward to derive the luminosity function of FRII sources 
from the model. Peacock (1985) and Dunlop & Peacock 

(1990) have used a large data base of deep radio obser- 



vations to determine the luminosity function of all radio 
galaxies and quasars. Their free modelling approach allows 
the shape of the luminosity function to be constrained but 
does not explain why it has this shape. We show that the 
model described here is in agreement with their results and 
allows some conclusions as to the origin of the form of the 
luminosity function. 

The median linear size of luminous sources in the LRL 
sample appears to be shorter than that of low luminosity 
sources (see Figure ^). Because of the correlation between 
luminosity and redshift in flux-limited samples this effect 
can also be interpreted in terms of a decrease in median lin- 
ear size, D me d, with redshift (e.g. Masson 1980 and Macklin 
1982). If interpreted in this way, this effect can be modelled 
with a power law, D me d oc (1 + z) nn . Several groups have 
tried to determine no and find very different values. Oort et 
al. (1987), Singal (1988), Kapahi (1989) and Subrahmanian 
& Swamp (1990) find n D ~ 3 while Eales (1985) and Neeser 
et al. (1995) derive no ~ 1-3. Neeser et al. (1995) showed 
that the strong evolution found in some samples could be 
caused by a selection effect in which FRI sources are mis- 
takingly included in samples of FRII sources leading to an 
overestimate of the median linear size at low redshift. How- 
ever, the observation that the maximum and the median lin- 
ear sizes of high redshift and therefore luminous sources are 
shorter than those of their lower redshift and less luminous 
counterparts must be explained by any model predicting the 
radio source distribution in the P-D plane. 

In Section 2 we construct samples of radio sources pre- 
dicted by various source distribution functions and compare 
their distribution in the P-D plane with that of the LRL 
sample. The luminosity function of FRII sources predicted 
by the best-fitting models is derived in Section 3. In Section 
4 we investigate the cosmological evolution of the median 
linear size of FRII objects as predicted by our model. A 
discussion of the implications of the model for the cosmo- 
logical evolution of the radio source progenitor population is 
presented in Section 5. Throughout this paper we are using 
Friedmann world models with A = and H = 50 km s _1 
Mpc" 1 . 



2 THEORETICAL SOURCE DISTRIBUTION 
IN THE P-D PLANE 

In this section we describe the construction of theoretical 
P-D diagrams and the method employed to compare them 
with the observed distribution of FRII sources in the P-D 
plane. 



2.1 The source distribution function 

The model presented in KA and KDA for the radio lumi- 
nosity and linear size evolution of a given radio source de- 
pends on a variety of source and environment parameters. 
If the density distribution of the gaseous environment of the 
source, p x , is approximated by a King (1972) profile, 



px = Po 



1 + 



(a o ) 



(1) 
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Figure 1. The source distribution of the LRL sample in the P-D plane. Left: Q = 1, right: Q = 0. Crosses: FRII-type radio galaxies 
and quasars, circles: FRI-typc radio galaxies. The binning of the P-D plane used to compare model predictions to the LRL sample is 
shown. 
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Figure 2. The source distribution of the LRL sample in the P-z plane. Left: Q = 1, right: Q = 0. Crosses: FRII-type radio galaxies 
and quasars, circles: FRI-typc radio galaxies. The binning of the P-z plane used to compare model predictions to the LRL sample is 
shown. 



where r is the radial distance from the centre of this distri- 
bution, these defining parameters are: The jet power, Q , 
the initial opening angle of the jet, 9, the central density of 
the density distribution of the source environment, p a , the 
core radius of this density distribution, a a , and the exponent 
of the power law describing this distribution, /3. Since it is 
not straightforward to measure 9 directly in observations, 
we will use the aspect ratio of the cocoon, R, i.e. the ratio 
of the length of the cocoon and its width, which is related 
to 9 by c 2 /9 2 w 4R 2 in the case of a ram-pressure confined 
cocoon, where c 2 is a known constant (see KDA) . 

We assume that both, the energy density of the par- 
ticles in the cocoon and that of the energy density of the 
magnetic field in the same region have relativistic equations 
of state which corresponds to case 1 in KDA. It was shown in 



that paper that the specific choice of the equation of state 
in both cases does not significantly influence the shape of 
the resulting evolutionary tracks of a given source through 
the P-D diagram. Note also, that our choice implies equipar- 
tition between the energy density of the particles and that 
of the magnetic field in the whole cocoon. We also assume 
that the power law describing the energy spectrum of the 
relativistic electrons in the cocoon of all model sources ex- 
tends from a Lorentz factor 7 = 1 to 7 — > 00 and that the 
exponent of the power law is p — 2.25, implying mildly rel- 
ativistic bulk velocities in the jets with a Lorentz factor of 
roughly 2 (Heavens & Drury 1988). 

For a given choice of source and environmental param- 
eters, together with the source age, the redshift of the 
source, z p , and the viewing angle between the jet axis and 
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the line of sight, a v , the model developed in KDA yields 
a point on the P-D plane. In order to construct a complete 
theoretical P-D diagram one ideally needs to know the distri- 
bution function of each of the parameters mentioned above. 
Furthermore, the distribution functions of various param- 
eters may not be independent, e.g. the sources with more 
powerful jets are also situated in higher density environ- 
ments. For simplicity we assume initially that the distribu- 
tion functions required in the model are independent. The 
number density of FRII sources in the universe as a function 
of jet and environmental parameters and of source age, ti, 
can then be written as 



d 6 p 



= (Hot iVi(log 10 Q ) d(log 10 Qo) N 2 (R) dR 
x Ar 3 (log 10 p a ) d(log 10 p ) N 4 (a )da 
x N 5 (f3)df3 JV 6 (*,)#,, 



(2) 



where ptot is the total number density of FRII sources in the 
universe and Ni, where i runs from 1 to 6, are the distribu- 
tion functions of the jet and environmental parameters. All 
Ni are normalised to yield 1 when integrated over the entire 
range of the respective parameter. 

In the chosen cosmology (A = 0) the progenitors of 
radio sources of age ti observed at redshift z p became active 
at a cosmological epoch corresponding to 



[Zp + 1) 5 

for Q = 1 and 

z p + 1 



3 Ho 



I, 



- 1, 



(3) 



(4) 



l-Hoti{z p + l) 

for Q = 0. Let the probability for a progenitor to 'switch 
on' as a function of redshift, the 'birth function' of FRII 
sources, be N 7 (z s )dz s . The distribution function of the ages 
of radio sources, ij, observed at redshift z p then is 



N 6 (ti)dti = N 7 (z s )^-dti 
dti 



= N' 7 (t u z p )- 

for Q = 1 and 
N 6 (ti)dti = N 7 (ti,z p )H 



H 



■dti 



z p + l 



1-Hoti {z p + 1) 



dti, 



(5) 



(6) 



for Q — 0. If not stated otherwise we assume the remaining 
distribution functions Ni(A), with i = 1. . . 5, to be uniform; 



An 



Ni (A) dA = f ( Amax > Amin) d, 1 . 
where 



1 5 Amin ^ A ^ Ajjiax 

; A < A m in or A max < A 



(7) 



(8) 



Ne{ti)dti in equations (^) and is the fraction of all 
FRII sources in the universe with an age ti at a cosmo- 
logical epoch corresponding to redshift z p . Of these only 
the sources located at coordinate distances between R r(z p ) 
and R r(z p + dzp) are observable. Assuming that the survey 
which we will use to compare our model predictions with 



subtends a solid angle of Q and that the highest redshift 
at which progenitors become active is z max , this fraction is 
found to be 



dV(z p ) 



Q [R r{z p )f d(R r) 
Q/3 [R r(z ma x)] 3 

3 (z max + l) 3 (V z p + 1 - 1)^ 



2 V (^ + !) J (Vw + l-l)' 
in the case of Q = 1 and 

dV{z p ) _ 3 {Zmax + l) 3 

V(z ma x) z 3 /l I I z ) 3 

v ' ^max \ ' 2 max I 

z% (l + ±z p ) 2 {l + z p + \z 2 p ) 



■ dz p , (9) 



dZr, 



(10) 



for Qo = 0. 

From the jet and environmental properties it is possible 
to determine the linear size of a given radio source for any 
source age, ti, and the observed linear size is determined 
by projecting onto the plane of the sky using a distribution 
function sin a v da v for viewing angles between a v and a v + 
da v . 

With these results the relative number of sources within 
the population of given jet and environmental parameters 
observed at redshift z p and at a viewing angle a v are de- 
termined. Computationally, this is achieved by dividing the 
range over which a property A varies into small intervals 
AA, which are then identified with the differentials dA in 
equation (jij). The model for the radio emission of the co- 
coon described in KDA us then used to convert the chosen 
source properties into a linear size and a radio luminosity 
at the frequency in the rest-frame of the source, P„( z 
corresponding to an observing frequency, v. The flux density 
which would be measured from the source is then 



S v 



(R ry(l + z p )- 



(11) 



where R r is the coordinate distance of the source. If this 
flux density is greater then the flux density limit of the ob- 
servational sample the source is included in the model sam- 
ple with the appropriate weight given by equation (^) and 
multiplied by dV (z p ) /V {z max ) sina„ da v . 

Objects in flux-limited samples of observed radio 
sources span a large range in redshift. To obtain radio lu- 
minosities for these objects at the same frequency in the 
rest-frame of the individual sources it is usually assumed 
that the spectrum of radio galaxies is well represented by a 
power law, Pi oc v~ a , where the spectral index, a, is deter- 
mined by measuring the flux of the radio galaxy at a second 
frequency. For the LRL sample this second frequency is 750 
MHz. With this assumption equation ( pd| ) can be inverted 
to give 



Pi = S v (R rf (1 + z p ) 



(12) 



In order to mimic this procedure we calculate the radio lumi- 
nosity of every source in the model samples at 178 MHz and 
750 MHz from the model of KDA and derive the 'observed' 
spectral index, a. With this and using equation ( |l^ ) we then 
calculate the 'observed' radio luminosity, P'(y), which is 
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then used to determine the location of the source on the P-D 
plane. Note, that Pi is independent of the value of Q while 
the flux density calculated from equation ([H]) depends on 
the chosen cosmology. The sets of sources within the model 
samples are therefore different for different cosmologies. 

Unless most of the radio emission from large sources 
is concentrated in the hot spots, it can be distributed over 
a large projected area on the sky. Even if the total flux of 
such sources is above the flux limit of a given sample, the 
sources may not be included in the sample because their 
surface brightness is below the detection limit of the tele- 
scope used to compile the survey from which the sample 
is drawn. In the model presented here we assume that the 
surface brightness of a given source is equal to its total lumi- 
nosity divided by the projected area of sky covered by the 
cocoons. If the calculated surface brightness is below the 
surface brightness limit of the 6C sample, which was used 
to check the completeness of the LRL sample (Riley 1989), 
the source is not included in the model sample. Using this 
selection criterion we may omit some sources from our model 
samples which should be included because their radio emis- 
sion is less smoothly distributed and their hot spots may 
be bright enough to be detected. However, in all the mod- 
els presented here, the 'number' of sources excluded because 
of insufficient surface brightness was of order 10 -4 at most, 
which is negligible when compared to the number of sources, 
10 -2 or higher, in the relevant bin which were included in 
the model sample. 



2.2 Model parameters 

The range over which the source parameters used in the 
model vary can in principle be deduced from observations, 
however, the constraints obtained in this way are incomplete 
and we will show later that there may also be considerable 
variation with redshift. 

For the jet powers of FRII sources with redshifts below 
8 = 1 Rawlings & Saunders (1991) find 10 37 W < Q rs < 
10 40 W from minimum energy arguments and Sl = 0. For 
Q, = l these limiting values of Q rs are somewhat lower 
but even for z = 1 the correction factor necessary is less 
than 2. KDA point out that another correction due to an 
underestimate of the expansion work done by the cocoon by 
Rawlings & Saunders (1991) and to the potential presence of 
thermal material in the cocoon is necessary. Their equation 
(18) assumes that the rate at which energy is lost by the 
cocoon as expansion work as predicted by the model of KA 
is equal to half the jet power found by Rawlings & Saunders 
(1991). Because of the method used by Rawlings & Saunders 
(1991) it is more accurate to assume the rate at which energy 
is stored in the cocoon as predicted by KA to be equal to 
Q rs /2. This implies that Q rs is an overestimate of the jet 
power. When the additional energy of the thermal particles 
required by the model in KDA is also taken into account, 
we find that the correction to Q rs is of order unity. We will 



than very luminous objects. We therefore introduce an ex- 
ponential distribution of jet powers, Q \ 



therefore initially use Q 0} . 



10 3 ' W and Q a . 



10* 



W. A further discussion of this point is deferred to Section 
3. 

The source distribution of the LRL sample in the P-z 
plane (see Figure ^|) implies that at every redshift there are 
more FRII sources of low radio luminosity in the universe 



-Wi(log 10 Q ) d (log 10 Q ) = 

\ e - A g ( lo glO Qo-log 10 < 







-A,(log 10 Q 



-l°glO Qa.min) 



d(log 10 Q o ), (13) 



where \ q varies from to 10 in steps of 1.0. For X q 







we use a uniform distribution of log 10 Q as described by 
equation (Q). 

For the aspect ratio of the cocoon, R, we use a uni- 
form distribution with 1.3 < R < 6.0, taken from Leahy 
& Williams (1984). The environments of FRII sources are 
either the IGM in clusters of galaxies or the atmospheres 
of the host galaxy itself in the case of small linear sizes 
or isolated objects. Canizares et al. (1987) find 4 • 10~ 23 
kg/m 3 < p < 5 ■ 10~ 21 kg/m 3 for the central density of 
individual galaxies. The core radius of the density profile 
in these objects is 0.01 kpc < a < 2 kpc with a median 
value d ~ 1.0 kpc. For clusters of galaxies Jones & Forman 
(1984) find 7 ■ 10" 25 kg/m 3 < p < 5 ■ 10" 23 kg/m 3 and 
30 kpc < a < 1100 kpc with d ~ 260 kpc. For simplic- 
ity we will calculate theoretical P-D diagrams for individ- 
ual sources and sources in clusters separately and we will 
also assume that the core radius, a , for all model sources 
is equal to the respective median value a . This implies 
N4(a ) da = 5(a — d ) da in equation (^). For the central 
density we use a uniform distribution with po.min = 5- 10 -23 
kg/m 3 and p . m ax = 5- 10 -21 kg/m 3 for isolated sources and 
po,min = 5 • 10~ 25 kg/m 3 and p ,max = 5 • 10" 23 kg/m 3 for 
sources in galaxy clusters. 

In some of the following models we allow for cosmolog- 
ical evolution of the central density, p . If the properties of 
the environments of the progenitors only depend on redshift 
but not on any other source property, i.e. the jet power, they 
should all evolve in very similar ways. The central density, 
po, will in this case be proportional to the density of the gas 
filling the universe at the time when the progenitor and its 
environment decoupled from the Hubble flow. In this sim- 
plified picture we therefore expect that p oc (1 + z) nd , with 
n d < 3. 

The model for the radio emission of the cocoons of FRII 
sources described in KDA requires an external density dis- 
tribution with p x oc r~P. For distances from the centre of 
the density distribution, r, greater than a few core radii, 
a , the approximation p x = po{r / a )~P provides a good 
fit to equation ([!]). However, within the core radius this is 
not the case and sources in the centre of galaxy clusters will 
spend a considerable part of their life time in this region. We 
therefore split the external density profile into three different 
regimes: p x = p a for < r < a D /2, p x = po/ (2r/a )" /3 ' /2 
for a /2 < r < 2a and p x = p (r/a ) ~^ for r > 2a . 
During its life time a given radio source is changing from 
one density regime into the next when the linear sizes of 
its cocoons, Lj, reach a /2 and 2a respectively. Let the 
cocoons of a given radio source reach a linear size of a /2 
each at an age of ti(Lj = a /2) = ti t o(Lj = a /2). Dur- 
ing this time the external density profile is assumed to be 
px = Po- If the source was located in a density regime with 
px = po (2 r/a )~' 5/ ' 2 while expanding to this linear size, 
its age would be ti^/ 2 {Lj = a /2) ^ ti t o(Lj = a /2). To 
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calculate the correct radio luminosity using the model pre- 
sented in KDA once the cocoons of the source become larger 
than a /2 and enter the second density regime, we have 
to use a ficticious source age ti^/ 2 (Lj > a /2) = ti(Lj > 
o /2) — tifi(Lj — a /2) + ti y p/ 2 (Lj = a /2) instead of the 
real source age ti(Lj > a a /2). In other words, once the linear 
size of the cocoons of the radio source become larger than 
a /2, we treat the source for the calculation of its radio lu- 
minosity as if it had been located in an environment with 
px = po (2 r/a„)~ fl ' 2 for a time equal to the fictious source 
age. Sources with cocoons larger than 2a entering the third 
density regime are treated in an analogous way. The aspect 
ratio of the cocoon, R, is assumed to remain constant when 
the source changes from one regime to the next. We assume 
a uniform distribution of /3 between 1.0 and 2.0. 

The viewing angle, a v , is distributed according to 
sina„ da v over the range to 7r/2 radians and for the source 
age we assume a maximum value of ti ymax — 10 9 years which 
is the upper limit of observed spectral ages (Alexander & 
Leahy 1987). For sources reaching the maximum age we as- 
sume that the jet ceases to supply the cocoon with energy. 
The subsequent drop in radio luminosity occurs fast because 
of the adiabatic expansion of the cocoons. Therefore we as- 
sume that the radio luminosity of a source of age ti trnax 
drops to zero instantaneously. 

The ratio of the number of progenitors becoming FRII 
sources with powerful jets to that of objects with weaker jets 
is given by equation ( |13[ ) and is constant for all redshifts. The 
distribution of observed sources within the LRL sample in 
the P-z plane (see Figure ^) therefore suggests that the total 
number density of FRII sources in the universe is higher at 
higher redshift. To allow for this effect we introduce a birth 
function of the form Nt(z s ) dz 3 oc (z s + l) n dz s . The distri- 
bution function of the age of sources observed at redshift z v 
then becomes 



Nr(ti,z p 



n + 1 



rrulx ~\~ 1) 1 



for 51 = 1 and 



(1 + %) 2 ~-^t, 



n + 1 

Zp + 1 

[1 - H ti (z p + 1) 



(14) 



(15) 



for fi = 0. In both cases we investigate values of n between 
and 10. Since there are no sources in the LRL sample 
beyond a redshift of roughly 2.1 we will use z m „ = 5 as 
the cut-off for the birth function. This ensures that even for 
Q — the first generation of radio sources which became 
active at a cosmological epoch corresponding to z 3 = 5 has 
reached the maximum source age at the cosmological epoch 
corresponding to z p = 2.1. 



2.3 x 2-ana ly s i s 

The comparison of the theoretical predictions for the distri- 
bution of radio sources in the P-D diagram with the LRL 



sample is done by splitting the P-D plane into 19 bins as 
shown in Figure oj. The number of FRII-type sources for the 
LRL sample in each of the bins is given in the first column 
of Tables |l| and |[ Note, that the FRI-type objects in the 
LRL sample are not included in the binning process. 

The way in which the P-D plane is binned and the num- 
ber of bins will influence how well particular model predic- 
tions agree with the observed distribution for the LRL sam- 
ple. Too many bins result in few sources in each occupied 
bin and it is then difficult within the x 2 analysis to dis- 
tinguish between occupied bins and those which are empty. 
As we have pointed out earlier, the maximum linear size of 
radio sources for each radio luminosity and the maximum 
radio luminosity within the whole sample are the strongest 
constraints on the cosmological evolution of the progenitor 
population. We have therefore chosen the bins in Figure [i] to 
delineate the source distribution of the LRL sample for large 
linear sizes and high radio luminosities without increasing 
significantly the number of bins. 

For each set of model parameters for the progenitor pop- 
ulation two 'samples' are calculated: one has external density 
profiles characteristic of isolated galaxies and the other den- 
sity profiles appropriate to cluster environments. In all cases 
the flux density and surface brightness limits match those 
of the LRL sample. These 'samples' are then normalised to 
contain in total the same number of FRII-type objects as 
the LRL sample and binned in the way shown in Figure |l| 
The normalisation determines the total number density of 
FRII sources in the universe, ptot, in equation (^|). 

The distribution of sources over the P-D diagram re- 
sulting from the model calculation is then compared with 
the one of the LRL sample using a \ test. The probability 
that the source distribution of the LRL sample is a repre- 
sentation of a universe with the chosen source distribution 
function, equation (^|), increases with decreasing \ 2 ■ 

Since the redshifts of all objects in the LRL sample are 
known it is also possible to compare the model predictions 
with the distribution of observed sources in the P-z plane. 
Figure |2| shows the objects in LRL in this plane and the 
17 bins used for the x 2 -test here. The binning of the P-z 
plane is complicated by the Malmquist bias. Most sources 
in the LRL sample and the model samples are close to the 
line defined by the flux limit and it is difficult to find bins 
in this plane that will quantify the scatter about this line 
properly. Ideally one would like to use the two dimensional 
Kolmogorov-Smirnov test to compare the models with the 
observed data in the P-z plane. However, Peacock (1983) 
showed that the two dimensional KS test is unreliable if the 
two independent variables, in this case P and z, are strongly 
correlated. Since the main interest in this analysis is the P-D 
diagram we will compare only the best-fitting models with 
LRL on the P-z plane using the much simpler \ 2 method. 



2.4 The cosmological evolution of the progenitor 
population 

24.1 Model A 

In the first model considered, model A, all distribution func- 
tions of source and environmental parameters except that of 
the jet power, Q , are assumed to be uniform as described 
by the 'top-hat' function, /, given in equation (fef). For the 
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A B C D E F G 



bin 


LRL 


gal. 


cl. 


gal. 


cl. 


gal 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


1 





0.73 


0.27 


0.66 


0.16 


0.54 


0.19 


0.61 


0.21 


0.61 


0.21 


0.34 


0.13 


0.35 


0.15 


2 


2 


6.78 


6.75 


5.95 


3.31 


4.81 


4.05 


6.06 


4.67 


5.05 


4.20 


3.26 


2.77 


3.51 


3.45 


3 


10 


16.46 


16.72 


13.24 


7.22 


10.59 


8.64 


11.22 


9.86 


8.96 


8.14 


7.23 


6.52 


8.09 


8.09 


4 


3 


2.41 


0.32 


1.44 


0.17 


0.73 


0.11 


0.33 


0.08 


0.60 


0.09 


0.46 


0.03 


0.35 


0.04 


5 


2 


0.67 


0.19 


1.11 


0.41 


1.00 


0.46 


1.28 


0.44 


1.24 


0.58 


0.90 


0.47 


1.35 


0.62 


6 


6 


5.06 


6.35 


6.72 


7.85 


5.55 


8.74 


8.24 


8.20 


7.23 


9.79 


7.32 


8.91 


9.46 


10.30 


7 


1(5 


12.64 


19.06 


12.26 


15.95 


9.17 


16.46 


11.49 


14.29 


11.04 


15.61 


12.65 


14.89 


13.81 


16.32 


8 





4.38 


1.21 


1.98 


0.49 


0.99 


0.28 


0.45 


0.13 


1.0 


0.18 


0.32 


0.00 


0.28 


0.00 


9 


5 


1.11 


0.32 


1.60 


0.91 


2.65 


1.38 


2.89 


1.34 


4.11 


2.35 


3.38 


2.07 


3.76 


2.12 


10 


23 


9.55 


9.22 


11.25 


14.52 


18.47 


21.71 


22.37 


20.50 


25.64 


27.25 


28.03 


31.20 


26.05 


26.60 


11 


32 


23.95 


29.86 


22.02 


28.89 


34.00 


39.41 


33.44 


35.57 


35.87 


32.39 


34.80 


34.43 


32.08 


26.19 


12 





11.03 


5.03 


3.28 


0.92 


2.10 


0.49 


0.73 


0.25 


1.33 


0.19 


0.12 


0.00 


0.25 


0.00 


13 


3 


2.74 


0.54 


2.88 


2.04 


5.04 


1.82 


4.24 


2.14 


5.87 


3.88 


4.16 


2.07 


4.52 


3.35 


11 


27 


19.45 


12.46 


19.79 


23.75 


21.51 


18.54 


21.48 


22.01 


20.86 


25.23 


25.94 


24.46 


22.47 


28.78 


15 


Xl 


12.39 


14.57 


11.79 


14.09 


24.83 


20.64 


18.30 


23.30 


12.80 


12.91 


13.55 


13.65 


16.32 


17.03 


16 





13.71 


20.24 


22.57 


17.74 


0.76 


0.09 


0.17 


0.06 


0.05 


0.02 


0.53 


1.37 


0.12 


0.00 


17 





0.00 


0.00 


4.58 


4.84 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


18 





0.01 


0.00 


0.12 


0.01 


0.25 


0.01 


0.18 


0.01 


0.20 


0.02 


0.01 


0.00 


0.25 


0.03 


19 





0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


x 2 
xl 




78.1 
78.3 


190.8 
171.3 


67.3 
66.3 


105.6 
58.2 


29.7 
23.0 


100.2 
26.5 


34.3 
12.8 


135.0 
24.9 


21.5 
12.0 


108.1 
11.4 


21.5 
7.7 


299.4 
17.1 


27.0 
6.9 


238.1 
12.7 



Table 1. The number of sources in each P-D bin as denned in Figure |lj for the LRL sample and models A to G for f2 Q = 1. gal. and cl. 
indicate the use of the isolated galaxy and the cluster density profile respectively, x 2 is the result of the x 2 -test for the respective model 
and x\ i s the result of the x 2 -test when omitting bin 4. 



A B C D E F G 



bin 


LRL 


gal. 


cl. 


gal 


cl. 


gal 


cl. 


gal. 


cl. 


gal. 




cl. 


gal. 


cl. 


gal. 


cl. 


1 





0.70 


0.15 


0.69 


0.17 


0.60 


0.19 


0.53 


0.21 


0.61 





.18 


0.39 


0.21 


0.37 


0.16 


2 


2 


6.23 


3.88 


5.98 


3.47 


5.17 


3.89 


5.35 


4.55 


5.13 


3 


.12 


3.73 


4.78 


3.72 


3.27 


3 


10 


13.01 


9.60 


11.25 


7.77 


9.77 


8.59 


10.17 


10.05 


9.16 


6 


.88 


8.41 


11.08 


8.34 


7.62 


4 


3 


1.14 


0.18 


0.73 


0.18 


0.64 


0.11 


0.27 


0.08 


0.54 





.13 


0.54 


0.04 


0.39 


0.05 


5 


2 


0.77 


0.15 


1.36 


0.42 


1.33 


0.40 


1.26 


0.39 


1.38 





.56 


1.03 


0.55 


1.62 


0.87 


6 


5 


6.09 


5.05 


7.86 


8.30 


7.75 


7.84 


7.65 


7.54 


7.81 


10 


.02 


8.55 


10.90 


10.37 


13.01 


7 


17 


13.84 


15.09 


13.55 


18.27 


13.51 


16.24 


12.07 


14.93 


12.41 


18 


.24 


14.40 


17.39 


15.37 


19.96 


8 





3.25 


0.98 


1.35 


0.62 


1.34 


0.31 


0.52 


0.17 


1.04 





.37 


0.36 


0.00 


0.31 


0.01 


9 


5 


1.50 


0.37 


2.36 


1.02 


3.06 


1.22 


3.03 


1.33 


3.88 


2 


.42 


3.17 


2.00 


4.23 


2.87 


10 


23 


13.92 


11.25 


15.88 


17.15 


20.78 


19.54 


23.03 


21.14 


24.11 


28 


.52 


28.08 


29.15 


26.72 


30.97 


11 


36 


33.08 


37.13 


28.69 


36.69 


37.56 


38.94 


38.97 


40.38 


34.38 


35 


.64 


34.77 


28.30 


30.14 


27.69 


12 





11.84 


7.30 


2.97 


1.45 


3.10 


0.82 


1.38 


0.46 


1.74 





.45 


0.16 


0.00 


0.27 


0.00 


13 


3 


1.88 


0.51 


3.70 


1.54 


3.49 


1.70 


2.61 


1.41 


5.25 


3 


.16 


3.08 


2.30 


4.79 


3.22 


14 


22 


14.55 


13.01 


19.87 


18.34 


16.50 


18.69 


16.58 


17.07 


20.98 


21 


.51 


22.58 


23.95 


22.48 


23.50 


15 


15 


19.08 


35.33 


22.97 


23.97 


18.18 


24.55 


19.49 


23.34 


14.36 


11 


.49 


13.76 


12.35 


13.65 


9.81 


16 





2.10 


3.07 


3.72 


3.75 


0.05 


0.02 


0.03 


0.01 


0.01 





.01 


0.00 


0.00 


0.00 


0.00 


17 





0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 





.00 


0.00 


0.00 


0.00 


0.00 


18 





0.01 


0.00 


0.15 


0.00 


0.16 


0.01 


0.12 


0.00 


0.20 





.01 


0.01 


0.00 


0.29 


0.03 


19 





0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 





.00 


0.00 


0.00 


0.00 


0.00 


x 2 




47.1 


180.4 


32.1 


80.9 


22.1 


100.4 


39.5 


139.6 


21.2 


75.9 


18.3 


227.0 


25.6 


205.9 


xl 




44.6 


138.3 


25.4 


38.3 


13.6 


28.5 


12.6 


27.5 


10.2 


14.4 


7.1 


18.1 


8.1 


17.6 



Table 2. The number of sources in each P-D bin as defined in Figure |l| for the LRL sample and models A to G for Ho = 0. gal. and cl. 
indicate the use of the isolated galaxy and the cluster density profile respectively. \ 2 is the result of the x 2 -test for the respective model 
and xl i s the result of the x 2 -test when omitting bin 4. 
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jet power the exponential distribution of equation ( |l3| ) is 
assumed. The birth function is assumed to be a power law 
of redshift and the distribution of source ages is therefore 
given by equations ( p^ ) and ([TH]) respectively. The best- 
fitting model parameters for both, the isolated galaxy and 
the cluster density profile of the progenitor environment are 
given in Table | for Q a = 1 and Table | for Q a = 0. The 
resulting relative source numbers in the P-D bins for the 
best-fitting model together with the value of the x 2_ test for 
this model are given in Table |l| for fl — 1 and Table |^ for 
tto = 0. 

For model A we note that the number of sources in 
the low luminosity bins are comparable to the LRL sample 
but in all cases there are too many powerful, large sources 
(bins 12 and 16). This implies that the median hot spot 
advance speed in the powerful sources is too high for the 
assumed set of distribution functions. The slope of the birth 
functions, equations ( |l4| ) and ([if]), is steep in all cases; the 
exponent n is in the range of 5 to 7. Although very luminous 
sources are therefore, in this model, located preferentially at 
a higher redshift than the less luminous ones, the increased 
energy losses of the relativistic electrons in the cocoon due to 
the higher energy density of the CMBR at high redshift do 
not sufficiently steepen the evolutionary tracks of the most 
powerful sources. 



2.4.2 Model B 

The hot spot advance speed of a given source depends cru- 
cially on the density distribution of the material the source 
is expanding into. An increase in the median external den- 
sity will therefore on average slow the linear expansion of the 
jets. The density of the progenitor environments at low red- 
shift, whether they are individual galaxy or cluster profiles, 
are constrained by observations. For model A most of the 
objects which have linear sizes larger than those of sources 
in LRL are luminous and therefore preferentially located at 
high redshifts. This trend strongly suggests that the density 
of the environments the progenitors are located in, increases 
with redshift, and therefore we have considered models in 
which upper and the lower limit of the range of the central 
density, p a , varies with redshift, according to 



Po,min — po^rairti^Zp — 0) i^Zp -f- 1) 
Po,max — Po^max^Zp — 0) (Zp ~\~ 1) 



(16) 



In these expressions p ,min(z p = 0) = 5 • 1CP 23 kg/m 3 
and po,max(z P = 0) = 5 ■ 10~ 21 kg/m 3 for galactic and 
po,min(zp = 0) = 5 • 1(T 25 kg/m 3 and p ,max(z p = 0) = 
5 ■ 10~ 23 kg/m 3 for cluster density profiles as before. The 
parameter rid is investigated in the range to 9. The best 
fit obtained for model B predicts a very strong evolution of 
p with redshift independent of the form of external density 
or Q . If the density of the progenitor environments is in- 
dependent of other source paramete rs a s we have assumed 
here, we expect < 3 (see Section 2.2). Model B is there- 
fore unphysical since it requires na > 9. 



2.4.3 Models C and D 

Although the exponential distribution of jet powers, equa- 
tion (|13|), ensures that there are many more progenitors de- 
veloping weak jets as compared to those producing strong 
jets, the results of model B suggest a persisting overestimate 
of the number of sources with very powerful jets in the model 
sample. We have therefore considered a variant of model B 
with a maximum jet power of log 10 (Q Oima:E /W) = 39.4 in 
the case of Sl = 1 and log 10 (Q Oima;E /W) = 39.7 for fl = 
instead of log 10 (Q o , mal /W) = 40. For the best-fitting model 
parameters the resulting model C includes almost no sources 
which are too luminous. Unless otherwise stated the follow- 
ing model samples are calculated with a maximum jet power 
of 10 39 ' 4 W for n o = 1 and 10 39 ' 7 W for 0, = 0. 

For the determination of the best-fitting set of model 
parameters in model C bin 4 was excluded from the x 2 -test 
since it alone contributes up to 76% of the total deviation 
of the model from the data of the LRL sample found with 
the x 2 ~test. The contribution of bin 4 is so large that it 
dominates the fitting procedure and any model for which 
the distribution of the LRL sample in the P-D plane may 
be a very good representation, could be rejected as a bad 
fit, if there is a strong deviation of LRL from the model 
sample in bin 4 alone. This suggests that some or all of the 
sources of the LRL sample found in bin 4 represent a class of 
objects distinct from the rest of the sample. These sources 
must be unusually old and/or located in very underdense 
environments. In the following we will exclude bin 4 from the 
X 2 -test when determining the best-fitting model parameters. 

Although model C represents an improvement over 
models A and B, particularly in the case of galactic den- 
sity profiles and Q = 0, the predicted redshift evolution of 
the central density, p , is again steeper than the physical 
limit in all cases. Another possible mechanism which could 
reduce the linear size of powerful jets at high redshift is the 
predominance of flat density profiles for the environments 
of progenitors at high z. This is consistent with the sce- 
nario of the cosmological evolution of the environments of 
the progenitors outlined above since we expect that at early 
stages of the evolution of regions which have decoupled from 
the Hubble flow, the density gradient within these regions 
is small. In model D we therefore introduce an exponen- 
tial distribution for j3 peaking at (5 — 1 which steepens for 
increasing redshift; 



N 5 ((3)d/3 

A/3 



A/3e -AgO?-/W) 
\p,o + (z P + 1)"" " 1 



dff 



(17) 



where A^,o and np vary from to 4 in steps of 1. Although 
model D improves the model fit significantly in the case of 
individual galaxy density profiles, the predicted evolution of 
the central density is still too steep with rid = 9. 



2.44 Model E 

For models A to D we have assumed that the maximum 
age, tmax, of all radio sources is 10 9 years. If t ma x decreases 
with increasing redshift, we expect the median linear size 
of radio sources to be smaller at higher redshift. Because 
of the strong correlation of redshift and radio luminosity in 
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D 





gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


(:/. 


gal. 


(:/. 


gal. 


cl. 


\ 


4 


4 


4 


3 


5 


4 


5 


5 


4 


4 


2 


1 


5 


5 


n 


U 


C 

o 


1 

1 


u 


f-~; 
u 


/I 
1 


I) 


•7 
i 


K 
U 


r. 
O 


Q 
O 


u 





I) 
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Q 


y 


Q 

y 
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u 


Q 
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1 


1 














71/3 














3 


3 














nt 


















1.0 


1.0 










n t , q 






















0.5 


0.6 






n q 


























2.0 


2.0 


xl 


78.1 


171.3 


66.3 


58.2 


23.0 


26.5 


12.8 


24.9 


12.0 


11.4 


7.7 


17.1 


6.9 


12.7 


(7 


0% 


0% 


0% 


0% 


15% 


7% 


75% 


10% 


80% 


83% 


97% 


45% 


99% 


75% 


-!ogio(P*ot) 


0.61 


1.87 


3.92 


5.31 


0.90 


2.88 


1.02 


2.42 


1.51 


2.17 


2.98 


5.18 


1.89 


1.71 



Table 3. The best-fitting model parameters for f2 = 1. To determine the best fir between a model and the source distribution of the 
LRL sample bin 4 in Figure |l| was omitted. The value of the x 2_ test for this case and the significance of the fit, a, are given for each 
model. The normalisation of the source distribution function, equation (0), is given as negative logarithm of ptot in units of Mpc - 3 . As 
before gal. and cl. indicate the use of the isolated galaxy and the cluster density profile respectively. 



D 





gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


gal. 


cl. 


A, 


4 


4 


4 


3 


4 


4 


5 


5 


4 


3 


2 


2 


5 


4 


n 


7 


7 


2 





3 


4 


7 


8 


1 


2 


2 


2 


5 


3 


n-d 






9 


9 


9 


9 


9 


9 


8 
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3 


3 














71/3 














1 


1 














nt 


















1.0 


1.5 










n t ,q 






















0.5 


0.6 






n q 


























2.0 


2.0 


xl 


44.6 


138.3 


25.4 


38.3 


13.6 


28.5 


12.6 


27.5 


10.2 


14.4 


7.1 


18.1 


8.1 


17.6 


a 


0% 


0% 


9% 


0% 


70% 


4% 


76% 


5% 


90% 


64% 


98% 


38% 


96% 


42% 


-1 °gio (Ptot) 


2.36 


3.19 


4.98 


6.51 


4.60 


4.76 


2.79 


2.66 


4.10 


4.06 


5.09 


5.50 


3.93 


5.44 



Table 4. The best-fitting model parameters for f2 = 0. To determine the best fir between a model and the source distribution of the 
LRL sample bin 4 in Figure |l| was omitted. The value of the x 2_ test for this case and the significance of the fit, a, are given for each 
model. The normalisation of the source distribution function, equation (^), is given as negative logarithm of ptot in units of Mpc -3 . As 
before gal. and cl. indicate the use of the isolated galaxy and the cluster density profile respectively. 



flux limited samples this may imply a decreasing median 
linear size with increasing luminosity. To investigate this we 
introduce 

tmax tmax.O 1) (-^) 

in model E, where t m ax,o is 10 9 years and nt varies from 0.0 
to 4.0 in steps of 0.5. Although the fit in model E is im- 
proved as compared to the previous models, the predicted 
evolution of the central density, p , is again very steep. The 
decrease in the maximum life time of sources at high red- 
shift which feature preferentially powerful jets, alone is not 
sufficient to reduce the median linear size of these objects as 
required by the observations. Note that, for the assumption 
of cluster density profiles, model E of all model samples cal- 
culated here, provides the best fit in both cosmologies with 
the observed LRL sample. 



2.4.5 Models F and G 

To make progress and obtain a model which fits the distri- 
bution of the LRL sample in the P-D plane without invok- 
ing an unphysical evolution of the density of the progenitor 



environment, we relax the assumption that all jet and envi- 
ronment parameter are independent of each other. 

The model of the evolution of the radio luminosity of 
powerful radio galaxies described by KDA assumes that the 
power of the jet, Q , is constant over the entire life time 
of the source. Together with the uniform maximum age for 
these objects assumed here, this implies that the total energy 
transported from the AGN to the large scale structure, Etot, 
during t max is linearly proportional to Q , which may not 
necessarily be the case in radio galaxies. If the powerful jets 
exhaust their energy supply faster than their weaker coun- 
terparts, e.g. tmax oc Q nt ' q , than this linear dependence is 
replaced by E to t oc Q a nt ' q and we expect the median linear 
size of the powerful sources to decrease. For model F we take 



tmax — tmax,0 



~ n t,q 



(19) 



where nt, q is varying from 0.0 to 0.9 in steps of 0.1. We also 
take the maximum jet power, Q ,max, to be 10 39 ' 7 W for 
Q — 1 and 10 40 W for Q = which improves the fit at the 
highest luminosities (bins 14 and 15). For the assumption 
of individual galaxy density profiles the model fit improves 
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to more than 97% significance for both cosmologies and the 
model predicts E to t oc Q o 5 . For cluster density profiles the 
fit is far less good. 

Rawlings & Saunders (1991) find that the power of the 
jet is roughly comparable to the Eddington luminosity of 
the black hole in the AGN powering the jet at z ~ 1. Kor- 
mendy & Richstone (1995) find a correlation of the mass of 
the central black hole and the mass of the bulge in nearby 
galaxies. Since we model the density profile of the progenitor 
environment with a power law (equation [ij), the total mass 
of the progenitor is proportional to the central density of 
this distribution, p . This may imply that the power of the 
jet in radio galaxies depends on the mass and the density of 
gas in the progenitor object. In model G we therefore take 



Po,min — Po,min('2'p — 0) 



po^raax — Po^maxi^Zp — 0) 



Qo 



Qo 



Qo,': 



(20) 



where p , mi „(2 p = 0) = 5 • 1CT 23 kg/m 3 and p 0>max {z p = 
0) = 5 ■ 10~ 21 kg/m 3 for galactic and p ,min(z p = 0) = 
5 • 1(T 25 kg/m 3 and p ,max{z p = 0) = 5 • 1(T 23 kg/m 3 for 
cluster density profiles as before. The maximum age of all 
sources is assumed to be 10 9 years and the maximum jet 
power is set to 10 39 4 W for n o = 1 and 10 39 ' 7 W for n o = 
as before. The exponent n q is varied from 0.0 to 4.5 in 
steps of 0.5. The significance of the fit for the assumption of 
galactic density profiles is comparable to that of model F, 
while for cluster density profiles the fit improves compared 
to the last model. 



2.4-6 Comparison of models 

Models F and G together with the assumption of galactic 
density profiles for the environment of the progenitors rep- 
resent the best fit with the source distribution of the LRL 
sample in the P-D plane in both cosmologies of all the model 
calculations presented here. For the assumption of cluster 
density profiles model E provides a better fit than either 
model F or G. However, the very steep evolution of the 
central density, p , with redshift predicted by model E is 
unphysical and we will therefore restrict attention in the 
following section to models F and G with the assumption 
of galactic density profiles for the environments of the pro- 
genitor objects. The fact that the best-fitting parameters for 
models F and G are very similar, if not completely identical, 
for the two cosmologies investigated is remarkable given the 
higher radio luminosity required of sources to be included in 
the model sample and the higher maximum jet power used 
in the case of Q a — 0. The source with the highest redshift in 
the LRL sample, 3C9, is located at roughly z = 2.1. The de- 
rived linear size and radio luminosity for an observed source 
at 2 = 2.1 with measured angular size and radio flux density 
are a factor 1.6 and 2.6 respectively greater for 0. o = than 
for Q. — 1. The resulting difference in the source distribu- 
tion of the LRL sample on the P-D plane can apparently be 
accounted for in model F and G simply by the inclusion of 
sources with higher jet power in the model sample for the 
case fi = 0, without changing the model parameters. 

To determine which of the two models is in better agree- 



ment with the LRL sample we now compare the source dis- 
tribution of the model samples with the best-fitting model 
parameters as given in Tables ^ and ^| in the P-z plane with 
that of the LRL sample. The results of the x 2 -test for the 
binning of the P-z plane as introduced in Figure ^| are sum- 
marised in Table ^. Here again we have omitted the three 
sources of the LRL sample in bin 4 of the binning of the 
P-D plane. Note here, that the low value for the maximum 
jet power used in models F and G implies effectively a max- 
imum radio luminosity for sources in the model samples. 
Together with the flux limit this means that the agreement 
of the model samples with the LRL sample in the P-z plane 
will be poor at the highest redshifts and also for the high- 
est luminosities. If we restrict attention to lower redshifts 
and lower radio luminosities by omitting bins 11, 12 and 13 
from the x 2 -test, we find good agreement between the model 
samples and LRL. See also the following section for further 
discussion of this point. In both cosmologies the agreement 
between model G and the LRL sample is better than that 
between model F and LRL. 



3 THE FRII LUMINOSITY FUNCTION 

Using the results found in the previous section for the source 
distribution function, equation (Q), it is now possible to de- 
rive the luminosity function of FRII sources, p(P v ,g). Fol- 
lowing the definition of Peacock (1985), p{P v , z) is the num- 
ber of radio sources per comoving volume measured in Mpc 3 
and per unit interval of log 10 (P„) with P v measured in W 
Hz -1 sr _1 . This luminosity function is identical to the source 
distribution function divided by log 10 (P„) but it is a func- 
tion of the radio luminosity, P u at frequency v and redshift, 
instead of the source and environmental parameters. To be 
able to compare our results with those of Dunlop & Peacock 
(1990) we use their observing frequency 2.7 GHz and calcu- 
late the 'observed' spectral index between 2.7 GHz and 1.4 
GHz. The luminosity function is then determined for 10 bins 
of equal width in log 10 P2.7 from 23.0 to 28.0 in the same way 
as the model samples in the previous section. Note, that we 
are now interested in all radio sources in the whole universe 
and not only in those which are observable to us. We there- 
fore do not have to multiply equation (g) by equation ^ or 
( p"o[ ) respectively. 

Figures ^ and ^ show the resulting luminosity function 
for models F and G with the assumption of galactic density 
profiles in the environments of the progenitors. The flatten- 
ing of the function for the lowest luminosities and the steep- 
ening at the high luminosity end are caused by the limits 
placed on the jet power, Q , and the density of the progen- 
itor environment. Sources of lower luminosity than approx- 



imately P2.7 = 2 ■ 10 24 W Hz -1 sr _1 are likely to have low 
jet powers which makes their jets susceptible to turbulence. 
Since FRI sources are not included in the model, the space 
density of low luminosity sources predicted here is a lower 
limit. 

Almost all of the sources in the LRL sample are close 
to the observational cut-off introduced by the flux limit of 
the sample (see Figure The one source in LRL in bin 11, 
the quasar 3C196, and the most distant source in LRL, the 
quasar 3C9, in bin 13 are clearly peculiar exceptions to this 
observation since they are far too luminous for their respec- 
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Table 5. The number of sources in each P-z bin as defined in Figure g for the LRL sample and models F and G for o = 1 on the left 
and Q = on the right, gal. and cl. indicate the use of the isolated galaxy and the cluster density profile respectively, x 2 i s the result 
of the complete x 2 -test for the respective model and x 2 is the result of the restricted x 2 -test which omits bins 11, 12 and 13. cr is the 
significance of the fit using the restricted x 2 -test. 
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Figure 3. The FRII luminosity function for ii 
and dotted: z=1.5. 
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1. Left: model F, right: model G. Solid line: z=0, dashed: z=0.5, dot-dashed: z=1.0 



tive redshifts. The remaining sources in the LRL sample in 
bins 12 and 13 are difficult to reconcile with the very low 
relative source numbers in the model samples for models F 
and G, particularly in the case of f2 = 1. 

By introducing rather low maximum jet powers in the 
previous section it was possible to reduce the number of 
sources of high luminosity in the model samples but this also 
leads to a lack of sources in the model samples with lumi- 
nosities and redshifts comparable to the highest luminosities 
and redshifts in the LRL sample. When included, the num- 
ber of sources with the highest jet powers at high redshift 
is overpredicted by the model. This implies that, although 
the simple power law assumed for the birth function of radio 



sources, equations (|14|) and (|l^) , provides good model fits at 
low redshift it is not a good fit to the 'true' birth function at 
high redshift. The 'true' birth function must flatten or even 
turn over at high redshift to explain the comparatively small 
number of sources at these cosmological epochs. This is sup- 
ported by a comparison of our results with those of Dunlop 
& Peacock (1990) presented in Figure |B|. The flattening and 
the turn-over of the luminosity function predicted by their 
free-modelling approach occurs at redshifts comparable to 
those at which our model is overpredicting the number of 
sources with powerful jets when these are included in the 
model samples. 

Dunlop & Peacock (1990) used an extensive data base 
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Figure 4. The FRII luminosity function for f2 = 0. Left: model F, right: model G. Solid line: z=0, dashed: z=0.5, dot-dashed: z=1.0 
and dotted: z=1.5. 
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Figure 5. The FRII luminosity function as a function of redshift. Left: P 2 .7 = 10 26 W Hz" 1 sr" 1 , right: P 2 .7 = W 27 W Hz -1 sr" 1 . 
Solid line: model F for Q = 1, dashed: model G for O d = 1, dot-dashed: model F for Q = and dotted: model G for f2 = 0. The error 
bars show the range of the luminosity function for Q = 1 as predicted by the free-modelling approach of Dunlop & Peacock (1990). 



of observed sources with a flux limit at their observing fre- 
quency at 2.7 GHz of two orders of magnitude lower than 
that of LRL at 178 MHz as used in the previous section to 
constrain and normalise the model. The agreement between 
their results and the prediction of our models, particularly 
model G, is therefore remarkable. 

The form of the cosmological evolution of the radio 
source population can be dominated either by space den- 
sity evolution, or by luminosity evolution. The evolution of 
the FRII population as predicted by models F and G is dom- 
inated by the space density evolution of radio sources. This 
is a result of our assumption of independence of source and 
environment parameters. The ratio of sources with high jet 
power and those with low jet power is the same at each red- 
shift and only the total number of radio sources is a function 
of z. Objects of similar jet power are also located in similar 
environments and have therefore similar radio luminosities 
regardless of their redshift. Exceptions to this rule are large 
sources which are already affected by the inverse Compton 
losses of the relativistic electrons in their cocoons. These 
sources have lower luminosities at higher redshifts because 
of the higher energy density of the CMBR at these epochs 



and show therefore some luminosity evolution. In the case 
of pure density evolution or pure luminosity evolution the 
slope of the luminosity function will stay constant at all red- 
shifts. Measuring the slope of the luminosity function pre- 
dicted by our models between P2.7 = 10 25 W Hz -1 sr" 1 and 
P2.7 = 10 27 W Hz" 1 sr" 1 in Figures | and | we find that this 
slope decreases from —1.2 at 2 = to —0.9 at 2 = 1.5 in the 
case of model F. For model G the slope decreases from —1.5 
at 2 = 0.0 to —1.2 at 2 = 1.5. This slight flattening of the 
luminosity function with increasing redshift was also found 
by Dunlop & Peacock (1990) when they used a combined 
density-luminosity evolution model to fit their data. 

Despite the luminosity evolution of large sources pre- 
dicted by our models the overall evolution of the FRII pop- 
ulation will be dominated by density evolution because of 
the set-up of our models. Fainter samples of FRII sources 
indicate that the evolution of the radio luminosity function 
for sources of lower luminosity is different from that of high 
luminosity sources (e.g. Dunlop & Peacock 1990). Some form 
of luminosity evolution therefore must take place. To model 
this behaviour of the radio luminosity function we would 
have to introduce a dependence of the birth function on the 
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jet power. This may also make it possible to avoid the low 
maximum jet power imposed on models F and G. However, 
this dependence would lead to an increase in the number of 
model parameters and since the LRL sample includes only 
the most luminous radio sources at any given redshift, the 
constraints on these additional parameters would be weak. 
We have therefore not investigated such models and our re- 
sults presented here apply mostly to the high luminosity end 
of the radio luminosity function. 



4 LINEAR SIZE EVOLUTION 

In order to determine the positions of a radio sources with 
given jet and environment parameters in the P-D plane we 
have already calculated the linear sizes of these sources. It 
is therefore straightforward to determine the median length 
of radio sources, D me d, as a function of redshift. 

If fitted by a power law, D me d oc (1 + z) nD , we find 
rin = —1.3 for model F and no = —1.1 for model G in 
both cosmologies. The model data is fitted well by a power 
law (correlation greater than 0.97) to z = 1.5. For higher 
redshifts the small number of sources in the model samples 
effectively prevents the determination of a meaningful value 
of D me d- Both values for ud are within the error range of the 
value found by Neeser et al. (1995) from a statistical analysis 
of the sources in the LRL and the 6C sample (e.g. Eales 
1985). There are no FRI-type objects in our model samples 
and our result can therefore be taken as a confirmation of 
the large selection effect pointed out by Neeser et al. (1995) 
in the derivation of the linear size evolution of radio sources. 

When considering all radio sources in the universe we 
find no = —0.4 for both models in both cosmologies. Also, 
the median linear size at z = is a factor 5 greater when all 
sources are taken into account as compared to the median 
linear size of the sources in the flux-limited model sample. 
This is consistent with larger sources being excluded from 
a flux limited sample because they are less luminous. The 
steeper linear size evolution in a flux-limited sample shows 
that this bias increases with increasing redshift. 



5 DISCUSSION 

In the models considered in this paper we have assumed the 
birth function, i.e. the fraction of the radio source popula- 
tion becoming active as a function of cosmological epoch, 
to increase monotonically with redshift. Since the central 
black hole powering AGNs must form at some cosmologi- 
cal epoch, one expects the birth function to turn over at 
some redshift and our model will predict too many sources 
at higher redshifts. This effect is indeed present in our mod- 
elling and led us to introduce a lower maximum jet power 
since only sources with the highest jet powers and therefore 
highest radio luminosities are observable beyond the turn- 
over redshift. We find that our model fits the observational 
data well out to a redshift of about 1.5. The turn-over in 
the birth function, and therefore the peak in the radio lu- 
minosity function, must therefore occur beyond, but not far 
from, 2 = 1.5. This is consistent with the results of Dunlop 

6 Peacock (1990). 

The FRII objects in the LRL sample represent the 



brightest radio sources in the universe at their respective 
cosmological epoch. They must therefore contain the most 
powerful AGNs and their cosmological evolution puts im- 
portant constraints on the formation of structure in the uni- 
verse. The median value of the jet power of sources in our 
model samples with the same flux limit as LRL as a function 
of redshift is well fitted by a power law and we find 
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for Q, = 0. The boundaries of the range of central densities 
of the progenitor environment, p , are constant for model F 
while they are proportional to for model G. We therefore 
expect the median value of the central density as a function 
of redshift in the model samples to be constant in model 
F and to increase with redshift as determined by equations 
@ and (§§) in model G. However, the introduction of a 
flux limit may result in a deviation from this and indeed we 
find 
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in the case of Q = 0. 

Equations (^) and imply that the total mass of the 
progenitors of the most luminous radio galaxies in model G 
increase quickly with increasing redshift. Assuming that the 
gas in these objects is distributed according to equation (|l]) 
with (3 = 2, the total mass of the gas in the progenitor out 
to a radius R x is given by 



Mtot = 7T p a ( — — arctan — \ 
\a a J 



(25) 



For the observed galactic density profiles, which we have 
used in this analysis, Canizares et al. (1987) use R x = 
50o o = 50 kpc for the largest radius at which X-ray emission 
is detected in their objects. In both cosmologies, we derive 
total masses for the progenitors within this radius of order 
10 9 M for model F and of order 10 10 M for model G at 
z = 0. Out to redshift 1.5 this mass increases insignificantly 
in model F for both cosmologies. However, in model G we 
find Mtot at z = 1.5 to be of order 10 14 Mq for Q = 1 and 
of order 10 15 M Q for fi D = 0. These derived mass concen- 
trations for the radio source progenitors within a radius of 
50 kpc are much higher than in any observed object. This 
may rule out model G. 

Rawlings & Saunders (1991) showed that the radiated 
luminosity of AGNs is roughly equal to the power of the jets 
associated with them. Since at least the radiated energy of 
AGNs is provided by accretion of material onto a central 
black hole, this implies that the jet power can be written as 



)„ = eL Bd = 1.25- 10 31 e^^W, 

Mp) 



(26) 



where Le<i is the Eddington luminosity of the central black 
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hole with mass Mbh- e is the fraction of the Eddington 
luminosity at which the AGN is supplying the jet with en- 
ergy. In order for accretion onto the central black hole to 
be possible, e < 1 is required. Theories of the formation of 
structure in the universe predict that the mass of the central 
black hole in any concentration of matter is proportional to 
the mass of the entire object, Mtot (e.g. Efstathiou & Rees 
1988). This is confirmed by the masses of central black holes 
derived from observations of nearby galaxies (Kormendy & 
Richstone 1995). With equation ( Eg ) this implies Q oc ep . 
Using equations (j5l]), (§), (§^), and (ji^) we find the fol- 
lowing expressions for e of the brightest radio galaxies as a 
function of redshift. 
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For increasing redshift e is increasing for model F while it 
is decreasing for model G. Rawlings & Saunders (1991) find 
the highest jet powers (~ 10 40 W) in sources at redshift 
1. If e — 1 in these sources, i.e. they are fuelling their jets 
at the Eddington limit, at z = 1.5, then the mass of their 
central black holes is of order 10 9 Mq. In this case we find 
e at z — for model F to be equal to 0.016 for o = 1 
and e = 0.006 if o = 0. For model G we find the rate 
at which energy is supplied to the jet to be highly super- 
Eddington at redshift (e = 85 for o — 1 and e = 162 for 
= 0). If e = 1 in the brightest radio galaxies at z — 
instead of at z = 1.5, then model G predicts AGNs with 
black holes in their centre which are sub-Eddington for all 
redshifts. However, this implies that the mass of the central 
black hole in the most powerful sources at z = 1.5 is equal to 
several 10 10 solar masses or higher. Although the existence 
of black holes with such large masses at high redshift can not 
be completely ruled out, the number of objects hosting such 
massive black holes is certainly very small (e.g. Efstathiou 
& Rees 1988). Note also, that in model F e oc Q while in 
model G e oc Qo 1 ■ This implies that the energy conversion 
mechanism driving the jets in radio galaxies is more efficient 
in sources with high jet powers then in those with weaker 
jets in model F while in model G exactly the opposite is 
true. 

In model F the life time of radio sources depends on 
the power of the jets. The median life time of the brightest 
radio sources is therefore also a function of redshift and we 
find tmax = 3.2 ■ 10 7 years for o — 1 and t max = 2.2 • 10 7 
years in the case of Q = at z = 1.5. This is consistent 
with observed spectral ages which usually do not exceed a 
few 10 7 years (Alexander & Leahy 1987). If the life time of 
the jets in radio galaxies is limited by the supply of fuel to 
their AGNs, then this implies that powerful jets are not only 
more efficient than weaker jets but also that they exhaust 
their fuel supply faster. 

The very high mass concentrations in the environment 
of the progenitors and the need for black holes of extremely 
high mass in the centre of the AGNs driving the most pow- 
erful jets in model G, suggest that this model is unphysical. 
However, KA point out that the general dynamics of a ra- 
dio source with given jet and environment parameters does 
not depend on p and a separately but only on the com- 
bination p a„ of these two quantities. The same is true for 



the radio luminosity of the source derived from the model of 
KDA. The only difference between the evolutionary tracks 
through the P-D diagram of two sources with the same value 
of poCLo but different core radii and central densities is caused 
by the different fractions of their life time they spend in the 
three density regimes introduced in Section 2.2. In the fol- 
lowing we will neglect the comparatively small effect these 
differences will have on the model parameters of model G 
and replace p a by p «o- With this we obtain p a a,o k Qi for 
model G. For the assumption that e = 1 for all values of Q 
we then find a oc Qo^ 3- ' 3 ' and p a oc Qq 6-3 ' 3 ^'' 3- ' 3 '. Using 
the redshift dependence of the jet power of the brightest ra- 
dio galaxies, equations ( pl| ) and (|5^), and assuming (3 = 2 
we note that p is now a constant while a oc (1 + z) 5,4 for 
o = 1 and a a oc (1 + z) 6 ' 2 for fi D = 0. Taking p a = 5 • 10" 22 
kg/m 3 , which is the median value of the distribution of p a 
for the galactic density profiles, and a — 1 kpc at z — 0, we 
find from equation (G5J) Mtot ~ 10 12 Mq within a radius of 
50 kpc at z = 1.5 for both cosmologies. This is comparable 
to the mass contained within the same radius in M87 found 
from globular cluster dynamics (Cohen & Ryzhov 1997) and 
from thermal X-ray emission (Nulsen & Bohringer 1995). 
Note, however, that in the case discussed here, the core ra- 
dius of the gas distribution of the environment, a a , increases 
to 141 kpc for o = 1 and 293 kpc for o — at redshift 
z = 1.5. Both values are much larger than those expected 
for the gas density profiles of individual galaxies and are 
more consistent with the gas distributions in galaxy clusters 
(Jones & Forman 1984). If e for the most luminous radio 
galaxies is not constant but decreasing with decreasing jet 
power, as is the case in model F, a will increase even more 
strongly with redshift while p may decrease with increasing 
redshift making the distribution of the gas in the environ- 
ments of radio galaxies at high redshift even more similar to 
that in galaxy clusters. 

Recently Best et al. (1998) showed that the host galax- 
ies of the most luminous radio galaxies at z ~ 1 are massive, 
highly evolved systems, presumably similar to the progeni- 
tors of the Brightest Cluster Galaxies (BCG). This is consis- 
tent with the tendency of the most luminous radio galaxies 
at high redshift to be located in richer environments than 
their low redshift counterparts (e.g. Hill & Lilly 1991). If we 
allow for the evolution of the core radius, a , of the density 
distribution in the environment of radio galaxies, our model 
G predicts a similar change in the environments of the most 
luminous radio galaxies from small, dense gas halos of indi- 
vidual galaxies at low redshift to more extended structures 
reminiscent of present day galaxy clusters at high redshift. 
This also explains the poor fit of our models with the LRL 
sample for the assumption of cluster-like density profiles in 
their environments at all redshifts. 

The strong evolution of the radio luminosity of the most 
luminous radio galaxies in the universe with cosmic epoch 
(see Figure ^ is explained differently in models F and G. 
For the assumption that the mass of the central black hole 
is proportional to the total mass of the radio source progeni- 
tor, model F predicts that the mass of the central black hole 
powering the most luminous radio sources is essential con- 
stant with redshift. The strong decrease in radio luminosity 
is therefore caused by a less and less efficient energy conver- 
sion process which expresses itself in the decrease of e with 
cosmic time. We have shown that model G is consistent with 
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a constant value of e and the decrease of the radio luminos- 
ity is then in this model caused by a decrease of the mass 
of the central black hole in the centre of the most luminous 
radio galaxies. In the latter scenario some additional pro- 
cess is needed to explain why the most massive black holes 
do not produce jets with high jet powers at low redshift. 
The faster virialisation of the material in large objects like 
galaxy clusters as opposed to smaller groups could prevent 
gas from reaching the centre of potential radio source hosts 
within these rich environments at low redshift and thereby 
depriving the black holes in these objects of fuel (i.e. Best et 
al. 1998, Ellingson et al. 1991). This may explain why radio 
sources in clusters at low redshift are usually of type FRI 
which implies that their jets are comparatively weak. 



with which the jets are fuelled. The alternative model G re- 
quires cosmological evolution not only of the central density 
but also of the core radius of the density distribution of the 
material surrounding radio galaxies. This is consistent with 
the observed change in the environments of radio galaxies 
with cosmic epoch. In this model, the decline of the radio 
luminosity of the most luminous sources is caused by the 
decreasing mass of the central black holes powering the jets. 
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6 CONCLUSIONS 

Based on the models for the evolution of the linear size and 
radio luminosity of powerful extragalactic radio sources of 
type FRII presented in KA and KDA, we have investigated 
the distributions of various jet and environment parameters 
and their evolution with redshift within the FRII source pop- 
ulation. The source distribution in the P-D plane predicted 
by our model is compared to that of the observed LRL sam- 
ple. We find that our model predicts an unphysically strong 
evolution of the gas density in the source environment with 
redshift, if we assume that all jet and environment parame- 
ters are independent of each other. The best fit of the model 
to the data in the P-D plane is achieved by assuming that 
the life time of radio sources or the shape of the density 
distribution of their environments depends on the power of 
their jets. 

Using this approach we find evidence that the giant 
sources in the LRL sample with linear sizes greater than 
1.5 Mpc, DA240, 3C236 and 3C326, constitute a class of ob- 
jects intrinsically different from the rest of the sample. They 
have to be extremely old and/or are located in extremely 
underdense environments. 

The luminosity function of FRII sources derived from 
the models is in good agreement with the results of Dunlop 
& Peacock (1990). We find evidence for a flattening of the 
luminosity function beyond z = 1.5. The evolution of the 
luminosity function is dominated by density evolution but 
pure density evolution is ruled out because of the inverse 
Compton scattering losses of the CMBR off the relativistic 
electrons in the cocoons of radio sources. This result only ap- 
plies to FRII sources with the highest radio luminosities at 
any given cosmological epoch which comprise the LRL sam- 
ple. For lower luminosities the cosmological evolution of the 
FRII source population is probably different and the models 
presented here have to be modified to allow for such a be- 
haviour of the radio luminosity function. Fainter complete 
samples of radio sources are also needed to better constrain 
the models at lower luminosities. The cosmological linear 
size evolution predicted by the models is consistent with the 
weak evolution derived by Neeser et al. (1995) from obser- 
vations. 

The best- fitting model predicting a correlation of the jet 
power with the source life time, model F, explains the decline 
of the radio luminosity of the most luminous radio galaxies 
with cosmic epoch in terms of a decrease in the efficiency 
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